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)  A  group  of  quadrature  formulae  for  end-point  singular  functions  is  presented  generalizing  classical 
end-point  corrected  trapezoidal  quadrature  rules.  The  actual  values  of  the  end-point  corrections  are 
obtained  for  each  singularity  as  a  solution  of  a  system  of  linear  algebraic  equations.  The  algorithm 
is  applicable  to  a  wide  class  of  monotonic  singularities  and  does  not  require  that  an  analytical 
expression  for  the  singularity  be  known;  only  the  knowledge  of  its  first  several  moments  and  the 
ability  to  evaluate  it  on  the  interval  of  integration  are  needed. 
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1.  Introduction 

When  a  numerical  quadrature  rule  for  a  class  of  singular  functions  is  derived,  it  is  usually 
a  generalization  of  some  technique  for  numerical  integration  of  non-singular  functions.  [3,  4, 
5]  use  Richardson-style  extrapolation  to  construct  quadrature  rules  for  singularities  of  the  form 
x1,  -1  <  -y  <  1.  [7,  8,  9]  generalize  the  Euler-Madaurin  expansion  for  integrands  with  certain 
kinds  of  singularities.  In  the  present  paper,  the  end-point  corrected  trapesoidal  quadrature  rules 
are  generalized.  While  the  standard  end-point  corrected  trapezoidal  rules  are  usually  derived  by 
means  of  the  Euler-Maclaurin  formula,  their  generalizations  for  singular  functions  are  obtained  in 
this  paper  as  solutions  of  certain  systems  of  linear  algebraic  equations  (see  Section  2).  Hie  rate 
of  convergence  of  resulting  quadrature  rules  is  investigated  analytically  for  a  fairly  broad  class  of 
monotonic  singularities  in  Sections  3-5.  However,  our  numerical  experiments  indicate  that  these 
quadratures  can  be  effectively  applied  to  a  considerably  broader  class  of  singularities  than  the  ones 
investigated  in  Sections  3-5  (see  Section  6). 

The  principal  advantage  of  the  approach  of  this  paper  is  its  generality  (see  Sections  3-6). 
Also,  since  almost  all  integration  weights  are  equal,  it  is  convenient  from  the  programming  point 
of  view  (see,  for  example,  [2]).  An  obvious  drawback  of  the  algorithm  is  the  fact  that  there  is  no 
simple  analytical  expression  for  the  corrections;  they  have  to  be  evaluated  numerically  for  each 
type  of  singularity,  and  for  each  subdivision  of  the  interval  of  integration,  unless  the  singularity  is 
either  of  the  form  log(x),  or  of  the  form  x'. 

2.  End-point  corrected  trapesoidal  rules  for  singular  functions 

A  right-end  point  corrected  trapezoidal  rule  T™  is  defined  by  the  formula 

W>  *=  k  •  (4f=i  +  2 /(*,))  +  *£>,■  /(*-/«)  (2-1) 

i«l  /-! 


where  /  :  (0, 1]  — *  f?1  is  an  integrable  function,  n,  m  are  a  pair  of  natural  numbers,  h  *=  £,  x;  =  for 
i  =  1, 2, ...,  n,  and  »  =  (oj,  0-2, ...,  am)  is  a  finite  sequence  of  real  numbers.  We  will  say  that  the  rule 
T”  is  of  right-end  order  k  >  1  iff  for  any  /  €  c*+1  [0, 1]  such  that  /( 0)  =  /'( 0)  =  ...  =  /(*)  (0)  =  0, 
there  exists  c  >  0  such  that 

I  T2(f)-£fW\<£-  (2.2) 

It  easily  follows  from  the  Euler-Maclaurin  formula  that  for  any  k  >  1,  there  exists  a  vector  a  = 
(ai,...,a*)  such  that  the  right-  end  order  of  the  formula  is  equal  to  k  (see,  for  example,  (6)). 
Suppose  now  that  the  function  /  :  (0, 1]  — *•  i?1  is  of  the  form 

f(x)  -  ■  s(x)  +  ip(x)  (2.3) 

with  4>.  €  c*[0, 1],  and  s  €  c(0, 1]  an  integrable  function  with  a  singularity  at  zero.  For  a  finite 

sequence  8  —  (81-82 ,...,/?m):  we  will  define  the  mapping  T"9  :  c(0.  l]  — *  I?1  by  the  formulae 


r:,!/)  =  r;(/)  +  /,^». 

i-l 


X.  = 


m 


(2.4) 

(2.5) 


t 
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We  will  use  expressions  T£g  with  appropriately  chosen  a,  0  as  quadrature  formulae  for  for  functions 
of  the  form  (2.3),  and  the  following  construction  provides  a  tool  for  finding  0  once  a  is  given. 

For  a  pair  of  natural  numbers  k,  n,  we  will  consider  the  following  system  of  linear  algebraic 
equations  with  respect  to  the  unknowns  0%,j  =  1,2,  ...,2k: 


2* 


y-i 


y.-i 

<K> 


(2.6) 


for  t  as  1,2,  ...,k  ,  and 


£  xf*"  •  *(Xj)  ■  07  *  f1  *(*)  •  ^*'1^  -  I?(e(x)  •  x*-1) 

•/o 


(2.6) 


for  i  **  k+  l,fc+2,  ...,2fc.  We  will  denote  the  matrix  of  the  system  (2.6),  (2.7)  by  XJ*,  its  right-hand 
side  by  Y’#w*  and  its  solution  (07,0%,  ...,0%t)  by  0%.  When  there  is  no  danger  of  confusion,  we  will 
omit  some  of  the  indices,  and  write  (0\,  02,.. .,02k)  and  0  or  0”,  instead  of  (0%,0%,...,0%k)  and  0% 
respectively.  The  use  ef  expressions  as  quadrature  formulae  for  functions  of  the  form  (2.3)  is 
based  on  the  following  theorem. 

Theorem  2.1. 


Suppose  that  the  systems  (2.6),  (2.7)  have  solutions  0" 
large  n,  and  that  the  sums 

B*nf 


(0%, 02,:; 02 it)  f°r  all  sufficiently 


(2.8) 


are  bounded  uniformly  with  respect  to  n.  Suppose  further  that  the  function  /  :  (0, 1)  — ♦  J?1  is 
defined  by  (2.3).  Then  there  exists  such  c  >  0  that 


I  T^n(f)-  ['  /(<)*!<  4- 

JO  n 


(2.9) 


for  all  sufficiently  large  n. 

Proof. 

Applying  the  Thylor  expansion  to  the  functions  we  obtain 

f(x)  =  P(f)(x)  +  Rk(<f>)(x)  •  s(x)  +  RkW, 


where 


i*0  i*0 

and  Rk(4>),RkW  are  such  functions  (0,1]  —  R1  that 


Rk(<l>)(  o)  =  (RkW)'(o)  = ...  =  (Rk(<(>)){k)(  o)  =  o, 

Rkm  o)  =  mm  o)  = ...  =  rnmHo)  *  o. 


(2.10) 


(2.11) 


(2.12) 

(2.13) 


Substituting  (2.10)  into  the  left-hand  side  of  (2.9),  we  obtain 

i  ir*(/)-  f l/m\< 

Jo 

-  f\a(x)  ■  R*(fa(x)  +  Rkmx))dx  | . 

Jo 

Due  to  (2.6),  (2.7), 

TS^P(f))  -  ['  P(f)(x)dx  -  0,  (2.15) 

Jo 

and  we  have  } 

I  T2r(f)-[  f(t)dt\< 

Jo 

I  TZ(a  ■  Rk(fa)  -  f  a(x)  ■  Rt(<l>)(x)dx  | 

Jo 

+  I  T2(Rk(rl>))  -  RtMWdx  1  (2-16) 

Jo 

+ 1  •  •ixj) 

)•  1 

+  W)(xy))  /?;i- 

Now  the  conclusion  of  the  theorem  follows  from  (2.2)  (2.8),  (2.12),  (2.13),  and  Lemma  A.2. 


3.  Convergence  rates  for  singularities  of  the  forms  x~*  and  log(x) 

For  the  remainder  of  this  paper,  </>\ ,  fa , ...,  fa  *  will  denote  functions  (0, 1]  —  R1  defined  by  the 
formulae 

<t>i{x)  =  x-1  for  i  -  1,2 (3.1) 
fa(x)  —  x,_1  ■  «(x)  for  i  =  k  +  1, k  +  2, ..., 2fc.  (3.2) 

Lemma  3.1. 

If  s(jc)  a®  x'7  with  -y  a  real  number  such  that  0  <|  *7  |<  1  then  the  functions  ^1,^2,  —,fak 
constitute  a  Chebychev  system  on  the  interval  (0, 1]  (i.e.  the  determinant  of  the  2k  x  2k  -matrix 
B  defined  by  the  formula  B.y  =  fa[tj)  is  non-zero)  for  any  2Ar  distinct  points  <1,(2,..., t2*  on  that 
interval. 

Proof. 

For  a  finite  sequence  Pi,P2,—,P2k  of  2k  real  numbers,  P  will  denote  a  mapping  (0,1]  —  B1 
defined  by  the  formula 

2k 

P(t)  =  (2.2) 

I*  1 

We  will  prove  the  lemma  by  showing  that  if  <i,t2,  —,t2k  are  distinct  points  on  the  interval  (0, 1]  and 
P(ti)  -  0  for  all » s»  1,2,..., 2 k  then  all  2k  coefficients  pi,P2,-,P2k  are  equal  to  zero.  By  repeated 


application  of  the  Rolle  theorem,  we  conclude  that  there  exist  k  distinct  points  pi,  pa, ...,  y*  on  the 
interval  (0, 1)  such  that  P^(y,)  *  0  for  all  i  «  1,2,  ...,k.  Differentiating  (3.2)  k  times,  we  obtain 

p{k)iVj)  =  p]~k  ■  •  fj"1  (3.3) 

i*i 

with  Wi  =  (*  -  1  +  Tf)(t  -  2  +  */)...(*  -  Hr  +  7)  ,  and  it  is  clear  that  w,-  ^  0  for  all  t  «  1,2,  ...,k .  Now 
it  follows  from  Lemma  A.3  that 

P»+i  *  Pm  *  —  *  Pik  *  0.  (3.4) 

Substituting  (3.4)  into  (3.2)  and  applying  Lemma  A.3  again,  we  see  that  p\  ■=  p2  *  ...  =  pk  m  0. 
Theorem  3.1. 

If  a(x)  =  X1  with  0  <|  n  |<  1  then  the  convergence  rate  of  the  quadrature  rule  T”B.  is  at  least 

It. 

Proof. 

It  immediately  follows  from  Lemma  3.1  that  the  matrix  of  the  system  (2.6),  (2.7)  is  non* 
singular.  We  rescale  the  system  (2.6),  (2.7)  by  multiplying  its  «  -  th  equation  by  (n  •  fc),_1  for 
i  —  1, 2, ...,  k,  and  by  (n  •  for  i  *  fc-f-  l.A-f-  2, ...,  2k,  obtaining  the  system  of  equations 


1°  ft  +  2°  ■  ft  +  ...  +  (2k)°  •  ftk 


0  an  _  (In  Ux  -  ron(l)) 


( r^W 


l1  ft  +  21  •  ft  +  ...  +  (2k)1  •  ftk  m  -  ^ 


l*-i  .  ^  +  2*-1  •  ft  +  ...  +  (2k)**1  •  1)} 

T  •  ft  +  2-1 -ft  +  ...  +  (2k)1  • 

1-r+i  .  +  21+1  •  ^  +  ...  +  (2k) 1+1  •  /?2"t  =  (Jo  xl+'dx 


(  f 1  T"i+k~ldr  —  Tnlr‘i+k~1)\ 

l1+*"i  •  /?;*  +  21+*-1  •  ft  +  ...  +  (2k)1+t-1  •  ftk  =  ^ ^ 

with  respect  to  the  unknowns  @ 2k'  and  we  will  denote  the  matrix  of  the  system  (3.5) 

by  J3*  ,  and  its  right-hand  side  by  Zk  .  Obviously,  Bt  is  independent  of  n,  and  \Zk  \  is  bounded 
uniformly  with  respect  to  n  due  to  Lemma  A.2.  Consequently,  the  solution  of  (3.5)  is  bounded 
uniformly  with  respect  to  n,  and  due  to  Theorem  2.1,  the  convergence  rate  of  T”fin  is  at  least  k. 

The  proofs  of  the  Lemma  3.1  and  Theorem  3.1  can  be  repeated  almost  verbatim  with  s(x)  — 
log(x)  instead  of  e{x)  =  X1,  resulting  in  the  following  theorem. 


If  a(x)  =  log(x)  then  the  convergence  rate  of  the  quadrature  rule  7^«  is  at  least  k. 

4.  Convergence  rates  for  more  general  singularities 
Theorem  4.1. 

Suppose  that  the  function  s  :  (0, 1]  — *  Rl  is  such  that 

xa'{x)  „ 

lim  — —  =  a  (4.1) 

*—o  a(x) 

with  0  <|  a  |<  1  .  Then  for  sufficiently  large  n,  the  matrix  A"*  of  the  system  (2.6),  (2.7)  is 
non-singular,  and  the  convergence  rate  of  T£g.  at  least  k. 

Proof. 

The  proof  below  effectively  consists  of  approximating  s(x)  on  the  interval  [xi>X2t]  by  the 
function  s  =  •  x°  and  using  the  Theorem  3.2. 

We  rescale  thesystem  (2.6),  (2.7)  by  multiplying  its  i-th  equation  by  (n-1)'*1  for «  «  1,2, 
and  by  (n  •  Jb)'"*"1  for  i  =  k  +  1,  k  +  2,  ...,2A*  ,  thus  obtaining  the  system 


1”  ■  0t  +  2°  ■«  +  ...  +  (2*)"  ■ «  - 

1-  ■  ff  +  2‘  •  +  ...+  (»)■  ■  J&.  $  *(*  ~  ^ W) 


1*-1  .  £J>  +  2*-1  -US  +  ...  +  (2*)*-1  •  = 

1°  •  4*4  '  +  2°  •  4*4  '  ^2+ 

®(Xi)  *(Xi) 


a-i  ~  (JbV-1^-!?^1)) 


(n  ■  A)**1 


,  a(X2t)  *„  _  (Jn  X1  ■a(x)dx~T”(x1  •  s(»))) 

•”  +  (  }  '  a(xi)  '8ik~  *(Xi M"*)1 


J*-1  .  fjxil .  +  2*-i  .  .  £"+ 

»(Xl)  »(Xi) 

.  ®(X2t)  _  (fn'x*->  s(x)dx-Ta"(x*-'-s(xm 

-  +  (  )  '  .(xi)  •(xt)-(»-*)*-1 


whose  matrix  we  will  denote  by  B£ .  Due  to  Lemma  A.2,  the  right-hand  side  of  (4.2)  is  bounded 
for  sufficiently  large  n,  and  by  applying  Lemma  A.5,  it  is  easy  to  see  that 

B£  =  B*  +  (4.3) 


with  (6g),k  =  1,2,...  a  sequence  of  2k  x  2k  matrices  such  that  lim,,™^  ||<5£||  =  0.  Now  it  follows 
from  Lemma  A.4  that  for  sufficiently  large  n,  ||(BJ)“1||  is  uniformly  bounded.  Thus,  the  solution 
of  (4.2)  is  uniformly  bounded  for  sufficiently  large  n,  and  due  to  Theorem  2.1,  the  convergence  rate 
of  is  at  least  k. 

Theorem  4.2. 

Suppose  that  k  >  1  and  a  function  s  :  (0, 1]  — ►  Z?1  is  such  that 

1..  s(x)  is  integrable  on  the  interval  (0, 1]. 

2..  S(x)  —  jjLj  is  integrable  on  the  interval  (0, 1]. 

3..  s  is  monotonic  in  some  neighborhood  of  0. 

4..  The  system  (2.6),  (2.7)  has  a  unique  solution  for  all  sufficiently  large  n. 

5..  The  convergence  rate  of  is  at  least  k. 

Then  the  convergence  rate  of  the  quadrature  rule  T^«  is  at  least  k  —  1. 

Proof. 


Again,  we  start  with  rescaling  the  equations  (2.6),  (2.7).  We  define  a  diagonal  2k  x  2k  -matrix 
D  by  the  formula 


(4.4) 


for  i  =  1,2,...,  2 k,  and  the  theorem  easily  follows  from  the  observation  that  A”*  =  A"*  •  D. 


Observation  4.1. 

The  combination  of  the  Theorems  3.1,  3.2,  4.1.  4.2  permits  one  to  apply  end-point  corrected 
trapezoidal  quadrature  rules  to  a  fairly  wide  variety  of  singular  functions,  including  linear  combi¬ 
nations  of  powers  of  x,  products  of  powers  of  x  with  log(x),  etc.  Further  generalizations  of  the 
quadrature  rules  T”gn  are  discussed  in  Section  7. 


5.  Asymptotic  behaviour  of  correction  coefficients  as  n  — *  oo 

An  obvious  drawback  of  the  expressions  as  practical  quadrature  rules  is  the  fact  that 
the  weights  •~,0%k  have  to  be  determined  for  each  value  of  n  by  solving  a  system  of  linear 

algebraic  equations  (albeit  a  small-scale  one),  in  this  section,  we  eliminate  this  problem  for  the 
cases  8(x)  =  x7  and  s(x)  =  log(x)  by  constructing  a  new  set  of  correction  weights  Ai,Aj,...,A 2k 
independent  of  n,  and  such  that  the  quadrature  rules  T£x  are  still  of  order  not  less  than  k. 

Lemma  5.1. 

Suppose  that  a  =  (01,02 . o m)  is  such  that  the  right-end  order  of  the  quadrature  formula 

T £  is  m.  Then  for  any  a  >  0  there  exists  c  >  0  such  that 

|  p‘+1  •  (T>(xa)  -  f  x°dx)  -  ?a+l  •  (T’(x-)  -  f  xadx)  |<  £  (5.1) 

Jo  Jo  P 


for  any  integer  p,  q  such  that  p  <  q  . 


iM '.  •«*  »It  *3  **  ,%«* , 


From  the  formula  (2.1),  it  is  easy  to  see  that 


72  (xa)  -  ['  x'dx  ■  j-  +  -^(-)‘  +  1  ^a(Ll.J.+  *)«  -  f'x'dx 

’  Jo  2q  qj^q*  qfz  q  Jo 

-<?)~-=£<=r+<i£.»(e=i±ir 


« '  *  2  *  qfz  «  Jo 


-  (f  )a+1 


with  iZp,  defined  by  the  formula 


Due  to  Euler-Maclaurin  formula, 


=  0((|) - (i)— ’)  -  0(;L. 


and  we  conclude  that 


J5<**)  '  £  *■*>  =  (2)*-1  ■  (!?(**)  -  »■<!*)  +  0(i  (5.5) 

which  is  equivalent  to  (5.1). 

Theorem  5.1. 

Suppose  that  lc,m  are  two  natural  numbers  such  that  k  <  m  -  1  and  that  a  =  (ai,ar2,...,am) 
is  such  that  the  quadrature  rule  7^  is  of  right-end  order  m.  Suppose  further  that  s(x)  *  i*'  with 
0  <|  7  |<  1  ,  and  that  the  coefficients  0i,0%,...,0%k  are  defined  as  solutions  of  the  system  (2.6), 
(2.7).  Then 

1..  There  exists  a  limit 

lim  0?  =  A,  (5.6) 

rt-^OC 

for  each  i  =  1,2,...,  2/:. 

2..  for  all  i  =  1.2,...,  2 k, 

Iff-*.  1=  O(^).  (5.7) 

3..  A,-  do  not  depend  on  m,  as  long  as  m  >  k  +  1. 

4..  The  quadrature  formulae  T£x  are  of  order  at  least  k. 


Suppose  that  p,q  are  two  natural  numbers,  and  p  <  q.  Obviously, 


0P  =  (St)*1  •  Zpk 

0"  =  (B*)-1  •  Z\  (5.8) 

P-F  =  B?(Zl-Zl). 

Due  to  Lemma  5.1,  there  exists  c  >  0  such  that 

113  -  311  <  <-^r).  (5-9) 

and  by  combining  (5.8)  and  (5.9),  we  see  that  for  some  s  >  0, 

<-k-  (5.10). 

Since  01,02,...  constitute  a  Cauchy  sequence,  taey  converge  to  some  limit  A  =  (Aj.  A2,...,  A2*), 
which  proves  1.,  and  2., 3., 4.,  easily  follow  from  (5.6), (5. 10), (5.1),  and  (5.5). 

The  proof  of  the  following  theorem  is  a  repetition  almost  verbatim  of  the  proofs  cf  the  Lemma 
5.1  and  Theorem  5.1. 

Theorem  5.2. 

If  under  the  conditions  of  the  Theorem  5.1,  the  singularity  »(*)  =  x"1,  is  replaced  with  the 
singularity  «(x)  *  log(x),  the  conclusions  1.  -  4.  of  the  Theorem  5.1  remain  correct 

Remark  5.1. 

For  singularities  of  the  forms  x"r,/og(x).  Theorems  5.1,  5.2  reduce  the  quadratures  to  a 
more  “conventionaT  form 

f1  f(x)dx  *  ro"A(/)  «  h  (£  /(*,.)  +  +  h  ■  £;  am/(x„-i+1)  +  h  jr  Xjfixj)  (5.8) 

Jo  -l  2  u  tz 

where  /  is  defined  by  (2.3),  h  *  £,x<  =  t  •  h  with  «'  «  1,2,  ...,n, \j  =  j fa  with  j  *=  1,2,..., 2 k,  and 

the  coefficients  A  j  are  independent  of  n. 

Remark  5.2. 

Theorems  3.1,  3.2,  4.1,  4.2,  5.1,  5.2  establish  that  under  certain  conditions  the  order  of 
the  quadrature  rules  T£gK.T£x  is  not  less  than  k.  where  2k  is  the  length  of  the  vectors  0n  » 
(fi1.*? . (A1.A2 . A2* ) .  Somewhat  more  detailed  versions  of  the  proofs  of  these  the¬ 

orems  show  that  the  order  of  these  quadratures  always  lies  between  k  and  I:  +  1.  and  is  in  some 
cases  equal  to  I:  +  1.  In  the  examples  presented  in  the  following  section,  we  note  the  cases  when 
the  latter  situation  occures. 

6.  Numerical  Results 

Algorithms  have  been  implemented  for  evaluating  the  quadratures  under  the  conditions 
of  the  Theorems  3.1,  4.1,  4.2,  and  the  quadratures  T£x  under  the  conditions  of  the  Theorems  5.1, 


5.2.  The  implementation  of  the  algorithm  for  evaluating  the  quadratures  T£g»  is  completely  stright- 
forward.  After  evaluating  the  matrix  A"k  and  the  right-hand  side  Y"k  of  the  system  (2.6).  (2.7), 
the  system  is  solved  by  means  of  a  LINPAK  Gaussian  elimination  subroutine,  and  the  resulting 
correction  weights  are  used  in  the  quadrature  formulae  (  see  formula  (2.4)). 

In  order  to  evaluate  the  coefficients  Aj,  A2,  ...A2*  for  singularities  of  the  forms  x'J,log(x),  we 
start  with  a  right-end  corrected  trapezoidal  rule  T"  of  order  20,  obtained  by  means  of  the  Euler- 
Maclaurin  formula.  Under  these  conditions, 

I  -  A.-  |<  O(^)  (6.1) 

for  all  i  =  l,2,...,2fc  (see  Theorems  5.1,  5.2),  and  for  reasonable  k,  the  convergence  of  to  A,-  is 
almost  instantaneous.  We  used  n  =  64  and  REAL  *16  arithmetic  to  obtain  15  significant  digits 
presented  in  the  examples  below. 

In  the  Tables  1*6,  the  coefficients  Ai,A2,...A2*  are  listed  for  singularities  of  the  forms 
log(x), xi, x~ 2, x», x~l ,x~i5,  and  k  =  2,3,4.  It  is  clear  from  these  tables  that  the  correction 
weights  Ai,  A2,  ...A2*  grow  rapidly  as  k  increases,  causing  increased  round-off  errors  and  making  the 
use  of  large  k  impractical  in  actial  computations. 

In  many  cases,  k  =  3  appears  to  be  a  reasonable  compromise,  resulting  in  roughly  4-th  order 
convergent  rules,  and  not  involving  unreasonably  large  correction  weights.  Combining  such  a  rule 
with  a  standard  4-th  order  end-point  corrections  at  the  right  end  of  the  interval  (see,  for  example, 
[l])  we  obtain  the  following  quadrature  formula 

f  f(x)dx  =  h  •  (£/(*,)  +  £ ’  /(*/)  £(/(**- 2)  “  4/(*»-i)  +  »/(*«)))  (6-2) 

®  I«1  jm  1 

with  h  =  A,  .  —  i  h  for  (  —  l,2,...,n,  and  \j  —  g  j  for  j  —  1,2.. ..,6.  For  the  six  singularities 
listed  above,  the  coefficients  Aj,  A2,  ...A6  can  be  found  in  the  middle  columns  of  the  Tables  1  -  6.  In 
the  Table  7,  convergence  results  are  presented  for  some  of  the  rules  T£g  and  T”x  with  k  =  3  and  T£ 
a  standard  4-th  order  end-point  correction.  Column  1  of  this  table  contains  the  numbers  of  nodes 
into  which  the  interval  [0, 1]  was  discretized.  In  all  cases,  the  integrand  was  of  the  form 

J{x)  =  (stn(23x)  +  cos(24x))  •  s(x)  +  (sin(2lx)  +  cos( 22x)),  (6.3) 

with  columns  2,  3,  4  containing  the  relative  errors  for  the  rule  T”x  with  «(x)  equal  to  log(x),x$, 
and  x~&  respectively.  Columns  5.  6  contain  relative  errors  for  the  rule  T"gn  with  s(x)  equal 

to  x~ »  •  log(x)  +  x~*  and  (log(x))2  respectively.  Finally,  column  7  contains  relative  errors  for 
the  standard  4-th  order  end-point  corrected  trapezoidal  rule  applied  to  the  function  (st'n(21x)  + 
cos(22x)),  presented  here  for  comparison. 

The  following  observations  can  be  made  from  Table  7.  and  are  typical  for  the  quadratures 

rpn  rpn 
1  a\'  ■*  Q0n  * 

1..  In  all  cases,  the  speed  of  convergence  is  roughly  the  same  as  for  the  end-point  corrected  trape¬ 
zoidal  rule  applied  to  the  function  (sin(21x)  +  cos(22x))  (column  7). 

2..  For  the  singularities  log(x)  and  xs ,  the  rules  T£x  display  a  typical  4-tb  order  convergence. 

ft 

3..  In  columns  3  -  5,  the  convergence  is  somewhat  erratic,  especially  for  a(x)  =  x“io,  in  which 
case  the  rule  seems  to  fail  after  n  =  320.  In  order  to  clarify  the  situation,  the  4-th  column  of 
the  Table  7  was  extended  with  n  =  2560,5120, 10240.20480.  The  extension  is  shown  in  Table 


8,  from  which  it  is  clear  that  the  bottom  part  of  the  column  4  in  Table  7  should  be  described 
as  erratic  convergence,  rather  than  a  failure  to  converge. 

4..  It  appears  from  the  5-th  column  of  Table  7  that  the  rule  T”gm  converges  quite  well  for  s(x)  ss 
( log(x ))2  ,  even  though  this  singularity  is  not  covered  by  any  of  the  theorems  of  the  present 
paper.  This  has  been  repeatedly  observed  when  the  singularity  s  is  monotonic,  and  in  such 
cases  a  version  of  Theorem  4.1  is  usually  fairly  easy  to  prove. 

7.  Generalizations  and  conclusions 

The  algorithm  of  the  present  paper  admits  several  strightforward  generalizations. 

1..  There  are  classes  of  singularities  not  covered  by  this  paper  for  which  some  versions  of  Theorems 
3.1  and/or  5.1  can  be  fairly  easily  proven.  The  convergence  of  the  quadrature  rule  for  one  of 
them  (  (j log(x ))2  )  is  demonstrated  in  the  preceeding  section. 

2..  Correction  nodes  xi,X2 ,...,X2*  do  not  have  to  be  equispaced.  A  different  distribution  of  nodes 
could  possibly  reduce  absolute  values  of  the  correction  weights  Xl5  Xj,  ...A2*  and  Pi,  fc,  ...fck, 
thus  improving  the  convergence.  However,  such  specialized  choice  of  correction  nodes  would 
have  to  be  performed  separately  for  each  singularity  s. 

3..  The  quadratures  can  be  easily  modified  to  handle  functions  of  the  form 

m 

f(x)  =  +  53  <i>i (ar)  ■  Si{x)  (7.1) 

•=i 

where  0, 4>\ , fa ,  <t>m  are  smooth  functions,  and  si,«2>  ...,sm  are  several  different  singularities. 
However,  it  is  easy  to  see  that  the  absolute  values  of  the  weights  /?"  tend  to  grow  very  rapidly 
as  m  increases,  and  the  author  doubts  the  practical  usefulness  of  such  rules  with  m>  2. 

4..  Quadrature  rules  T£x,  T”gn  have  fairly  obvious  analogues  in  two  and  three  dimensions.  How¬ 
ever,  the  proofs  of  multi-dimensional  versions  of  Theorems  3.1.  3.2,  4.1,  4.2,  5.1  are  somewhat 
more  involved  than  those  of  their  one-dimensional  counterparts.  These  results  will  be  reported 
at  a  later  date. 
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Correction  Weights  for  the  Singularity  s  =  log(x) 

k=2  k=3  k=4 

0. 16012984 1535717E+01  0.222876601846009E+01  0.309348340177712E+01 

-.33825585219 1949E+0 1  -.123121207006261E+02  -.310178837674079E+02 

0.3627 888464434 13E+Q 1  0.315796599730867E+02  0.136205915590327E-I-03 

-.134662835787181E+01  -.384039159001043E+02  -.314747480872421E+03 

0.226735045911525E+02  0.421505412761263E+03 

-. 526589398 196889E+01  -.328785403878733E+03 

0. 138801 167137067E+03 
-  .245552 1037 18723E+02 


Table  2 

Correction  Weights  for  the  Singularity  s  =  x$ 
k=3  k=4 

0.140373389574362E+01  0.176138469558481E+01 

-  .61062697 546597 1E+01  -.138211834485298E+02 

0. 14582 14411 89467E+02  0.545915011781337E+02 

- .  16396 1527863182E+02  - .  1 1735748454987 1E+03 

0.895228275571686E+01  0.150779019932162E+03 

-.  193573822942934E+01  - . 1 14778491 157932E-I-03 

0.476230959836121E+02 
-.829784263315958E+01 


Table  3 

Correction  Weights  for  the  Singularity  s  = 
k=3  k=4 

0.515686238420012E+01  0.788957615797699E+01 

-.371802541857231E+02  -.10 148391026933 1E+03 

0. 105027804788073E-I-03  0.498205235333950E+03 

-.138820513091881E+03  -.124177860454341E+04 

0.879943 1664 14350E+02  0. 1 75 10939935 8045E+ 04 

-  .2 16782 16536 1041E+02  -. 1419085 15209795E+04 

0.617986326801910E+03 
- . 1 12327 464963600E+03 


k=2 

0.333895462377735E+01 
-.103691855551396E+02 
0. 12388 1739056139E+02 
-.485794297425160E+01 


k=2 

0. 107622636973350E+01 
-.147247373014210E+01 
0.138293501775035E+01 
-.486687657341754E+00 


] 


Table  4 

Correction  Weights  for  the  Singularity  s  =  zt 


It— 2 

0. 120275440902998E+0 1 
- .  191590588964576E+01 
0.18902 152 1886824E+0 1 
-.677063738252464E+00 


k=3 

0.159679209113902E+01 
-.749191954692780E+01 
0. 1826 10902534004E+ 02 
-.210496743900247E+02 
0. 1 179979575 18642E+02 
-.261608415945115E+01 


k=4 

0.207566347869413E+01 
-.17681 779330607 1E+02 
0.723376732327415E+02 
- .  159272316539246E+03 
0.207188787576306E+03 
- .  15867 1438914327E+03 
0.660717800591290E+02 
-.11 5483695626909 E-f 02 


Table  5 

Correction  Weights  for  the  Singularity  s 


k=2 

0.245072941794592E+01 
-  .671933 1 1 1554788E+01 
0.775314064392466E+01 
-.298453894632270E+01 


k=3 

0.363 848721531386E+01 
-.239520591801675E+02 
0.653732223792864E+02 
-  .84042 1842099929E+02 
0.520805019262272E+02 
- . 12597968 130667 1E+02 


k=4 

0.538493619088162E+01 

-  .63685958687 4 130E+02 
0.30097935367452 1E+03 

-  .73 1 1686983 10276E+03 
0. 10 126 16552452 19E+04 
-. 80942528393743 1E+03 
0.548533467508443E+03 
-.627343688909117E+02 


Table  6 

Correction  Weights  for  the  Singularity  s  =  x~"& 

k=2 


0. 184790352839992E+02 
-.765330172344763E+02 
0.992955952836217E+02 
-.407416133331446E+02 


k=3 

0.321040181373839E+02 
-.289780736142061E+03 
0.895220650 1555 13E+03 
- .  12632677 1498802E+04 
0.841976333390815E+03 
-  .2 15752550553635E+03 


k— 4 

0.535441666101165E+02 
-  .846379842042659E+03 
0.456845469197409E+04 
- . 1212964 1603 1938E+05 
0. 179000593914225E+05 
- .  1 50266053430178E+05 
0.673612868399159E+04 
-.125506014574396E+04 


Examples  of  Convergence  of  Quadrature  Rules 


log{g) 

X* 

x~& 

-2 

X  * 

■log(x) 

+x~* 

(M*))2 

.503E-01 

.326E-01 

.157E+00 

.870E-01 

.105E+00 

.246E+00 

.412E-03 

.478E-03 

.189E-01 

.121E+00 

.590E-02 

.177E-01 

.546E-04 

.498E-04 

.964E-03 

.148E-02 

.219E-03 

.115E-02 

.102E-04 

.506E-05 

.422E-04 

.584E-04 

.651E-05 

.722E-04 

.992E-06 

.359E-06 

.146E-05 

.243E-05 

.434E-07 

.453E-05 

.803E-07 

.233E-07 

.359E-08 

.679E-07 

.180E-07 

.283E-06 

.600E-08 

•147E-08 

.958E-08 

.159E-08 

.241E-08 

.177E-07 

.425E-09 

.91  IE- 10 

•157E-08 

.616E-09 

•234E-09 

.109E-08 

Table  8 

Extended  Convergence  Results  for  s(x)  =  x~'& 


10 

20 

40 

80 

160 

320 

.157E+00 

.189E-01 

.964E-03 

.422E-04 

.146E-05 

.359E-08 

640 

1280 

2560 

5120 

10240 

20480 

.958E-08 

.157E-08 

.218E-9 

.269E-10 

.319E-11 

.334E-12 
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